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An effective algorithmic method is presented for finding the local conservation laws for 
partial differential equations with any number of independent and dependent variables. 
The method does not require the use or existence of a variational principle and reduces 
qq . the calculation of conservation laws to solving a system of linear determining equations 

similar to that for finding symmetries. An explicit construction formula is derived which 
yields a conservation law for each solution of the determining system. In the first of two 
papers (Part I), examples of nonlinear wave equations are used to exhibit the method. 
Classification results for conservation laws of these equations are obtained. In a second 
paper (Part II), a general treatment of the method is given. 

1 Introduction 



In the study of differential equations, conservation laws have many significant uses, par- 
ticularly with regard to integrability and linearization, constants of motion, analysis of 
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solutions, and numerical solution methods. Consequently, an important problem is how 
to calculate all of the conservation laws for given differential equations. 

For a differential equation with a variational principle, Noether's theorem jl2], [|, ||, 
^, [p|] gives a formula for obtaining the local conservation laws by use of symmetries of 
the action. One usually attempts to find these symmetries by noting that any symmetry 
of the action leaves invariant the extremals of the action and hence gives rise to a sym- 
metry of the differential equation. However, all symmetries of a differential equation do 
not necessarily arise from symmetries of the action when there is a variational principle. 
For example, if a differential equation is scaling invariant, then the action is often not 
invariant. Indeed it is often computationally awkward to determine the symmetries of 
the action and carry out the calculation with the formula to obtain a conservation law. 
Moreover, in general a differential equation need not have a variational principle even 
allowing for a change of variables. Therefore, it is more effective to seek a direct, algo- 
rithmic method without involving an action principle to find the conservation laws of a 
given differential equation. 

In Ref. (2| we presented an algorithmic approach replacing Noether's theorem so as to 
allow one to obtain all local conservation laws for any differential equation whether or 
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not it has a variational principle. Details of this approach for the situation of ordinary 
differential equations (ODEs) are given in Ref. 0. Here we concentrate on the situation 
of partial differential equations (PDEs). 

In the case of a PDE with a variational principle, the approach shows how to use the 
symmetries of the PDE to directly construct the conservation laws. The symmetries of 
a PDE satisfy a linear determining equation, for which there is a standard algorithmic 
method || [l4| to seek all solutions. There is also an invariance condition, involving just 
the PDE and its symmetries, which is necessary and sufficient for a symmetry of a PDE 
with a variational principle to correspond to a symmetry of the action. The invariance 
condition can be checked by an algorithmic calculation and, in addition, leads to a direct 
construction formula for a conservation law in terms of the symmetry and the PDE. This 
approach makes no use of the variational principle for the PDE. 

In the case of a PDE without a variational principle, the approach involves replacing 
symmetries by adjoint symmetries of the PDE. The adjoint symmetries satisfy a linear 
determining equation that is the adjoint of the determining equation for symmetries. 
Geometrically, symmetries of a PDE describe motions on the solution space of the PDE. 
Adjoint symmetries in general do not have such an interpretation. The invariance condi- 
tion on symmetries is replaced by an adjoint invariance condition on adjoint symmetries 
and there is a corresponding direct construction formula for obtaining the conservation 
laws in terms of the adjoint symmetries and the PDE. The adjoint invariance condition 
is a necessary and sufficient determining condition for an adjoint symmetry to yield a 
conservation law. 

In general for any PDE, with or without a variational principle, the approach of Ref. Q 
for finding all local conservation laws gave the following step-by-step method: 

(1) Find the adjoint symmetries of the given PDE. 

(2) Check the adjoint invariance condition on the adjoint symmetries. 

(3) For each adjoint symmetry satisfying the adjoint invariance condition use the direct 
construction formula to obtain a conservation law. 

If the adjoint symmetry determining equation is the same as the symmetry determining 
equation then adjoint symmetries are symmetries. In this case the given PDE can be 
shown to have a variational principle. Conversely, if a variational principle exists, then 
the symmetry determining equation of the PDE can be shown to be self-adjoint, so that 
symmetries are adjoint symmetries. 

In order to solve the determining equation for adjoint symmetries, one works on the 
solution space of the given PDE. On the other hand, in order to check the adjoint 
invariance condition, one must move off the solution space by replacing the dependent 
variable(s) of the given PDE by functions with arbitrary dependence on the independent 
variables of the given PDE. The same situation arises when checking for invariance of 
the action in Noether's theorem. 

These steps of the method are an algorithmic version of the standard treatment pre- 
sented in Ref. [T^ ] for finding PDE conservation laws in terms of multipliers. In particular, 
multipliers can be characterized as adjoint symmetries that satisfy the adjoint invariance 
condition and thus can be calculated by the step- by-step algorithm (1), (2), (3). 

In this paper (Part I) and a sequel (Part II) , we significantly improve the effectiveness of 
this method by replacing the adjoint invariance condition by extra determining equations 
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which allow one to work entirely on the solution space of the given PDE. Consequently, 
by augmenting the adjoint symmetry determining equation by these extra determining 
equations we obtain a linear determining system for finding only those adjoint symmetries 
that are multipliers yielding conservation laws. At first sight it is natural to proceed as in 
Ref. [Q by solving the adjoint symmetry determining equation and then checking which 
of the solutions satisfy the extra determining equations. On the other hand, all of the 
determining equations are on an equal footing, and hence there is no requirement to 
solve the adjoint symmetry determining equation first. Indeed, as illustrated later in the 
examples in this paper, it is much more effective to start with the extra determining 
equations before considering the adjoint symmetry determining equation. This is true 
even in the case when the given PDE has a variational principle. 

In solving the determining system one works completely on the solution space of the 
given PDE. Hence one can use the same algorithmic procedures as for solving symmetry 
determining equations in order to solve the conservation law determining system. In 
particular, existing symbolic manipulation programs that calculate symmetries can 
be readily adapted to calculate solutions of the conservation law determining system. 
Moreover, for each solution one can directly obtain the resulting conservation law by 
evaluating the construction formula working entirely on the solution space of the given 
PDE. 

The conservation law determining system together with the conservation law construc- 
tion formula give a general, direct, computational method for finding the local conser- 
vation laws of given PDEs. We refer to this as the direct conservation law method. As 
emphasized above, its effectiveness stems from allowing the calculation of conservation 
laws to be carried out algorithmically up to any given order by solving a linear determin- 
ing system without moving off the solution space of the PDEs. Compared to the stan- 
dard treatment of PDE conservation laws, the determining system solves a long-standing 
question of how one can delineate necessary and sufficient determining equations to find 
multipliers by working entirely on the solution space of the given PDE. Most importantly, 
by mingling the adjoint symmetry equations with the extra equations in the determining 
system, one can gain a significant computational advantage over the standard methods 
for finding multipliers. 

In Sec. ||we illustrate the direct conservation law method through classifying conserva- 
tion laws for three PDE examples: a generalized Korteweg-de Vries equation, a nonlinear 
wave-speed equation, and a class of nonlinear Klein-Gordon equations. The classification 
results obtained are new in that they establish the completeness of certain families of 
conservation laws which are of interest for these PDEs. These examples show how to cal- 
culate all conservation laws up to a given order and also how to determine which PDEs 
in a specified class admit conservation laws of a given type. 

In the second paper (Part II) , we present a general derivation of the conservation law 
determining system and construction formula, and we also give a summary of the general 
method. 
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2 Examples of conservation law classifications 

Here we illustrate the use of our direct conservation law method on three PDE examples. 
For each example we derive the conservation law determining system and use it to obtain 
a classification result for conservation laws. 

The first example is a generalized Korteweg-de Vries equation in physical form, for 
which there is no direct variational principle. The ordinary Korteweg-de Vries (KdV) 
equation is well-known to have local conservation laws of every even order |pjj| , which 
can be understood to arise from a recursion operator E3]. Using our conservation law 
determining system, we derive a direct, complete classification for all conservation laws 
up to second order for the generalized KdV equation. This example illustrates a gen- 
eral approach for finding and classifying conservation laws for non-variational evolution 
equations. 

The second example is a scalar wave equation with non-constant wave speed depending 
on the wave amplitude. This wave equation has a variational principle and admits local 
conservation laws for energy and momentum arising by Noether's theorem from time- 
and space- invariance of its corresponding action. By applying our conservation law 
determining system, we classify all wave speeds for which there are extra conservation 
laws of first-order and obtain the resulting conserved quantities. This example illustrates 
a general approach for classifying nonlinear evolution PDEs that admit extra conservation 
laws. 

The third example is a general class of nonlinear Klein-Gordon equations. The class 
includes the sine-Gordon equation, Liouville equation, and Tzetzeica equation, which 
are known to be integrable equations jlj with local conservation laws up to arbitrarily 
high orders, starting at first order. Through our conservation law method we give a 
classification of nonlinear Klein-Gordon equations admitting at least one second-order 
conservation law. As a by-product we obtain an integrability characterization for some 
of the equations in this class. This example shows a general approach to classifying 
integrable PDEs by means of conservation laws. 



2.1 Generalized Korteweg-de Vries equations 

Consider the generalized Korteweg-de Vries equation 

G = u t + u n u x + u xxx = (2.1) 

with parameter n > 0. This is a first order evolution PDE which has no variational prin- 
ciple directly in terms of u and which reduces for n — 1 , 2 to the ordinary KdV equation 
and modified KdV equation, respectively. Its symmetries with infinitesimal generator 
Xu = T) [& |l4| satisfy the determining equation 

= D t Tj + u n D x rj + nu n - 1 u x r ] + D 3 x tj when G = (2.2) 

where D t = d t + u t d u + u tx d Ul + u tt d Ut H and D x — d x + u x d u + ;/, , (',, + u tx d Ut H 

are total derivative operators with respect to t and x. The adjoint of Eq. (2.2) is given 

by 

= -D t u - u n D x uo - D 3 x uj when G = (2.3) 
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which is the determining equation for the adjoint symmetries lu of the generalized KdV 
equation. (Note, since these determining equations are not self-adjoint, the only solution 
common to Eqs. (2.2) and (2.3) is uj = rj = 0.) 

The generalized KdV equation (2.1) itself has the form of a local conservation law 

D t {u)+D x (-±-jU n+1 +u xx ) = 0. (2.4) 



We now consider, more generally, local conservation laws 

D t $* + DJ> X = (2.5) 



on all solutions u(t,x) of Eq. (|2.l|) . Clearly, we are free without loss of generality to 
eliminate any dependence on u t (and differential consequences) in the conserved densities 
$*, <5> x . All nontrivial conserved densities in this form can be constructed from multipliers 

A on the generalized KdV equation, analogous to integrating factors, where A depends 
only on t, x, u, and x derivatives of u. In particular, by moving off the generalized KdV 
solution space, we have 

D t $* + D x $ x = (u t + u n u x + u xxx )A Q + D x {u t + u n u x + u xxx )A 1 + ■■■ (2.6) 

for some expressions A Q , A-p . . . with no dependence on u t and differential consequences. 
This yields (after integration by parts) the multiplier 

+D X ($ X -T) = (u t +u n u x +u xxx )A, A = A ~D X A 1 + --- (2.7) 

where T = when u is restricted to be a generalized KdV solution. We now derive 
an augmented adjoint symmetry determining system which completely characterizes all 
multipliers A. 

The definition for multipliers A(t, x, u, d x u, . . . , d^u) is that (ut + u n u x + u xxx )A must 
be a divergence expression for all functions u(t,x) (not just generalized KdV solutions). 
This determining condition is expressed by 

= E u ((u t + u n u x + u xxx )A) 
= -D t A - u n D x A - D 3 X A + (u t + u n u x + u xxx )A u - D x ((u t + u n u x + u xxx )A dx J 

+ ■■■ + {-lfD p x ((u t + u n u x + u xxx )A 9S J (2.8) 

where E u = d u —D t d u —Dd u +D f Dd u +D 2 d u +■■■ is the standard Euler operator 
which annihilates divergence expressions. Equation (2.8) is linear in ut, 1H X , u txx , . . . , and 
thus the coefficients of Ut and x derivatives of Ut up to order p give rise to a split system 
of determining equations for A. The system is found to consist of the adjoint symmetry 
determining equation on A, 

= -V t A - u n D x A - D 3 x A (2.9) 
and extra determining equations on A, 

k=l 



o = (i - (-mA d « a 

0=(l-(-lf)A^. 
Here V t = d t - {uu x + u xxx )d u 
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k=q+l H K HJ 



q=l,...,p-l 



(2.10) 

(uu x + u xxx ) x d Un . + ■ ■ ■ is the total derivative operator 
which expresses t derivatives of u through the generalized KdV equation ( [2.l|) . Conse- 
quently, one is able to work on the space of generalized KdV solutions u{t, x) in order 
to solve the determining system (2^) and ( 2.10 ) to find A(t, x, u, d x u, . . . , d^u). The de- 
termining system solutions are the multipliers that yield all nontrivial generalized KdV 
conservation laws. 

The explicit relation between multipliers A and conserved densities <&' , <f> x for gener- 
alized KdV conservation laws is summarized as follows. Given a conserved density $*, 

we find that a direct calculation of the Ut terms in + D & x yields the multiplier 

equation 



{u t + 



DT 



(2.11) 



where E u = d u —D x d u +D x d u +• ■ • is a restricted Euler operator and T is proportional 
to ut + u n u x + u xxx and differential consequences. Thus we obtain the multiplier 



A = £?„(**). 



(2.12) 



Conversely, given a multiplier A, we can invert the relation (2.12) by a standard method 
[jl4| using Eq. (2.11) to obtain the conserved density 



dX uA(t, x, Xu, Xd x u, Xd x u, . . .). 



(2.13) 



From Eqs. ( [2.12] ) and ( 2.13 ) it is natural to define the order of a generalized KdV con- 
servation law as the order of the highest x derivative of u in its multiplier. Thus, we see 
that all nontrivial generalized KdV conservation laws up to order p are determined by 
multipliers of order p which are obtained as solutions of the augmented system of adjoint 



symmetry determining equations (2.9) and (2.1C). 

Through the determining system (2.9) and (2.10) we now derive a complete classifica- 



tion of all conservation laws (2.5) up to second order, corresponding to multipliers of the 
form 



A(t 

The classification results are summarized by the following theorem. 



(2.14) 



Theorem 2.1.1 The generalized KdV equation (2A_ ) for all n > admits the multipliers 

1 



A = 1, A = u, A = u x 



,71+1 



n+1 



(2.15) 



The only additional admitted multipliers of the form (2.1^) are given by 



A = tu — x, if n = 1, 



(2.16) 
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This classifies all nontrivial conservation laws up to second order for the generalized KdV 
equation for any n > 0. 

The conserved densities for these conservation laws are easily obtained using the con- 
struction formula ( 2.1 3| ) as follows. For the multipliers ( [2.15| ) we find, respectively, 



$ = u 
2 



,71 + 2 



(n + l)(?i + 2) 



DO, 



(2.18) 
(2.19) 
(2.20) 



where D x 9 is a trivial conserved density. In physical terms, if we regard u as a wave 
amplitude as in the ordinary (n = 1) KdV equation, then these conserved densities 
represent mass, momentum, and energy fig]. 



For the additional multipliers (2.1£) and (2.17) we find 
1 

-i 

2 



$* = — iw 2 — xu, if n = 1, 



1 



1 

12' 



G 



1^(9, if n = 2. 



(2.21) 
(2.22) 



Proof of Theorem 2.1.1 : For multipliers of the form ( 2.14 ) the adjoint symmetry 
equation ( |2.9| ) becomes 

= -D 3 X A - u n D x A -A t + A u (u n u x + u xxx ) + A u Ju n u xx + nu 71 ' 1 ^ 2 + u xxxx ) 



+A u ^(u n u xxx + 3nu n 1 u x u xx + n(n - l)u n 2 u x i + u xxxxx ) 
and the extra equations ( [2.10 ) reduce to 



= A 



(2.23) 

(2.24) 
(2.25) 



Note that Eq. (2.24) is a differential consequence of Eq. (2.25). 

We start from equation (2.25). Its term with highest order derivatives is u xxx A u u 
and hence A„, „, 



0. This yields that A is linear in u xx , 
A = a(t, x, u, u x )u xx + b(t, x, u, u x ). 



(2.26) 



Then the remaining terms in Eq. ( 2.25 ), after some cancellations, are of first order 

= b Ux - a u u x - a x . (2.27) 



We now turn to the adjoint symmetry equation (2.23) and separate it into highest 
derivative terms in descending order. The highest order terms cancel. The second highest 
order terms involve u XXXXl and these yield 



= Da = a x + a u u x + a u ^u xx . 



(2.28) 
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Clearly, Eq. ( 2.28) ) separates into 

a Ux = a u = a x = 0. (2.29) 

Hence b Ux = from Eq. ( gTgg ). 

The next highest order terms remaining in the adjoint symmetry equation ( 2.23] ) in- 
volve u xx . These terms yield 

= -a t ~ 3b xu + 3(nu n ~ 1 a - b uu )u x 

which separates into b uu = nu n ~ 1 a and a t = — '3b xu . Hence we have 



— |—aw n+1 + (6i0) - -xat)u + b 2 {t, x). 
n + 1 3 



(2.30) 
(2.31) 



This simplifies the remaining terms in Eq. (2.23) 




1 



b 2xxx - b 2 t - (bit - T, xa tt)u + b 2x u n + (- 

3 3 n + 1 



a t u 



n+l 



Clearly, the terms here can be separated, yielding 

(n - 2)a t = 0, 



2x x x i 



b 2 t = ~b 
bit - ^xa tt = -b^vJ 1 ^ 1 . 



(2.32) 

(2.33) 
(2.34) 

(2.35) 



From Eq. (2.35) there are two cases to consider. 
Case (i): n = 1 
Here Eq. yields 

a t = (2.36) 

and hence Eq. ( |2.35 ) becomes bit = ~b 2x . Since bi depends only on t, using Eq. (2.34) 
we find b 2xx — b 2t — 0, and thus 

b\ = b^t + 64, b 2 = —b^x + &5, 63 = const, 64 = const, 65 = const. (2.37) 

(2.38) 

(2.39) 



Finally, from Eqs. (|2.29|) and (|2.36| ) it follows that 

a = const. 

Therefore, through Eq. ( |2~2l| ), Eq. ( |2~3l| ), Eqs. ( p7| ) and ( ^38| ), we have 

A = -^au 2 + au xx + 63 — 652; + (64 + b$t)u, if n = 1, 

which is a linear combination of the multipliers ( |2.15| ) and ( 2.16| ) shown in Theorem 2.1.1 
Case (ii): n > 1 



In this case Eq. ( 2.35 ) clearly splits into 

bit = a tt = 0, 
b 2x = 0. 



Hence from Eq. fl2.4l| ) and Eq. ( |2.34| ) we have 

62 = const, 



(2.40) 
(2.41) 

(2.42) 
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and from Eq. ( 2.40 ) and Eq. Q2.29 ) we have 

a = ait + a 2l ai = const, a 2 — const, 
b\ = const. 



Therefore Eq. (2.31) becomes 



b = — — -{ait + a 2 )u n+1 + (bi - \xai)u + b 2 
n + 1 3 



(2.43) 
(2.44) 



(2.45) 



Finally, consider Eq. ( |2.33| ). If n 2 then we have at — and hence ai = 0. Otherwise, 
if n = 2 then we have no restriction on a\. Consequently, through Eq. (2.26), Eqs. (2.43) 
and (2.45), we find 



A = (ait + a 2 )u xx + 



,n+l 



(ait + a 2 ) + (bi - -xai)u + b 2 , if n > 1, 



(2.46) 



where oi = for all n > 2. This yields a linear combination of the multipliers (2.15) 
and Q2.17[) shown in Theorem |2.1.1[ □ 



2.2 Nonlinear wave-speed equation 

Consider the wave equation 

G = u tt -c(u)(c(u)u x ) x ={) (2.47) 

with a non-constant wave speed c(u). This is a scalar second-order evolution PDE, which 
has a variational principle given by the physically-motivated action 

S = J \{-ut 2 + c(ufu x 2 )dtdx. (2.48) 

Symmetries of the wave equation ( |2.47| ) with infinitesimal generator Xu = rj j^, Q satisfy 
the determining equation 

= D 2 t n - c(u) 2 D 2 x r) - 2c(u)d ' (u)u x D x T) - 2c(u)c' (u)u xx r] 

~(c(u)c"(u)+c'(u) 2 )u x 2 r] when G = 0. (2.49) 

Since Eq. (2.47) is variational, the determining equation ( 2.49| ) is self-adjoint and hence 
the adjoint symmetries of the wave equation are the solutions of Eq. ( [2.49] ). 

For any wave speed c(u), the wave equation (2.47) clearly admits time and space 
translation symmetries 

n = u t ,n = u x , (2.50) 

which are symmetries of the action. In particular, the action is invariant up to a boundary- 
term 

XS = J ^D t (-u t 2 + c(u) 2 u x 2 )dtdx (2.51) 
under time-translation Xu = u t , and 

XS= J ^D x (~u t 2 + c(u) 2 u x 2 )dtdx (2.52) 

under space-translation Xu = u x . By Noether's theorem, combining the invariance of the 
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action with the general variational identity 



XS* = 



((utt - c 2 u xx - cc'u x 2 )rj + D t (-u t rj) + D x (c 2 u x ri)jdtdx, Xu = i], (2.53) 



we obtain corresponding local conservation laws D t $ + D X ^> X = on all solutions u(t, x) 
of the wave equation (2.47). The conserved densities <& x are given by 

*' = \ut 2 + \c{u) 2 u x 2 , = -c(u) 2 u x u t (2.54) 
for the time-translation, and 

$* = u x u t , = ~^u t 2 - \c{u) 2 u 2 (2.55) 

for the space-translation. These yield conservation of energy and momentum, respectively. 

We now consider classifying wave speeds c(u) that lead to additional local conservation 
laws of first-order for the wave equation (2.47). (Through Noether's theorem all first-order 
conservation laws correspond to invariance of the action 5* under contact symmetries.) 

For any local conservation laws 

D t $* + D x $ x = (2.56) 



on all solutions u(t,x) of the wave equation (|2.47| ), we are clearly free without loss of 
generality to eliminate any dependence on utt (and differential consequences) in the con- 
served densities $ , <& x . All nontrivial conserved densities in this form can be constructed 

from multipliers A on the wave equation, analogous to integrating factors, where A has 
no dependence on uu and its differential consequences. In particular, by moving off the 
wave equation solution space, we have 

D t <S> 1 + D x $ x = (u« - c(u)(c(u)u x ) x )A + D x (u tt - c(u)(c(u)u ;r ) a: )A 1 + • ■ • (2.57) 

for some expressions A , A 1; . . . with no dependence on w tt and differential consequences. 
This yields (through integration by parts) the multiplier 

D t & +D X ($ X -T) = (u tt -c(u)(c(u)u x ) x )A, A = A - D X K X + • • ■ (2.58) 

where Y — when u is restricted to be a wave equation solution. We now derive an 
augmented symmetry determining system which completely characterizes all multipliers 
A. 

Multipliers A are defined by the condition that (« lf — c{u)(c{u)u x ) x )k is a divergence 
expression for all functions u(t,x) (not just wave equation solutions). We restrict at- 
tention to A of first-order, depending on t,x,u,Ut,u x . This leads to the necessary and 
sufficient determining condition 

= E u ((u tt - c(u)(c(u)u x ) x )K) 

= D 2 A - c 2 D 2 x A - 2cc'u x D x A - 2cc'u xx A - (cc" + c' 2 )u x 2 A 

+ {u tt - c 2 u xx - cc'u x 2 )A u - D t ((u tt - c 2 u xx - cc'u x 2 )A Ut ) 

-D x {(u tt - c 2 u xx - cc'u x 2 )A u ) (2.59) 
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where E u = d u 



D t d u 



D x du 



~D t D x d Uf 



is the standard Euler operator 



which annihilates divergence expressions. Equation (2.5£) is quadratic in u u and linear 
in uut and Uttx, and thus splits into separate equations. We organize the splitting by 
considering terms in G 2 , G, D X G, D t G, and remaining terms with no dependence on uu 

and differential consequences. The coefficients of D X G, D t G, and G 2 are found to vanish 

as a result of A(i, x, u, ut, u x ) being first-order. The other coefficients in the splitting do 
not vanish. 

This leads to a split system of two determining equations for A(t, x, u, ut,u x ), consisting 

of 



= V t A - c{u) 2 D x A - 2c(u)c'(u)u x D x A - 2c(u)c' (u)u xx A 



(c{u)c"(u) 



which is the symmetry determining equation (2.49) on A, and 

= 2A„ 



V t A 



D *K 



c'(u) 2 )u x 2 A 
(2.60) 

(2.61) 



which is an extra determining equation on A. Here T> t — dt + u t d u + u tx d Um + (c(u) 2 u xx + 
c(u)c' (u)u x 2 )d Ut +- ■ ■ is the total derivative operator which expresses t derivatives through 
the wave equation (2.47). Consequently, one is able to work on the space of wave equa- 
tion solutions u(t,x) in order to solve the determining system ( 2.60| ) and (2.61) to find 
A(t, x, u, ut, u x ). The determining system solutions are the multipliers that yield all non- 
trivial first-order conservation laws of the wave equation ( 2.47 ). 

The first determining equation ( 2.6C ) shows that all multipliers are symmetries of the 
wave equation, while the second determining equation ( 2.61 ) provides the necessary and 
sufficient condition for a symmetry to leave the action ( [2.48 ) invariant up to a boundary 
term. This is a consequence of the one-to-one correspondence between symmetries of the 
action and multipliers for nontrivial conservation laws of the wave equation (2.47) as 
shown by Noether's theorem fTi) . 

The explicit relation between multipliers A and conserved densities $ , $ x for first- 
order conservation laws (2.56) is summarized as follows. Given a conserved density <£"', 

we find that a direct calculation of the u tt terms in D t & + D x $ x yields the multiplier 
equation 



D& 1 



(u tt - c(u)(c(u)u 2; ) :c )£ , „ t ($ t 



DV 



(2.62) 



where E Ut = d Uf is the truncation of a restricted Euler operator, and T is proportional 
to u t t — c(u)(c(u)u x ) x and differential consequences. Thus we obtain the multiplier 



A = E Ut (d> t ). 



(2.63) 



Conversely, given a multiplier A, we can invert the relation ( [2.63 ) using Eq. ( 2.62 ) to 
obtain the conserved density 



& = 



where 



d\((ut~ut)A[\u + (l-\)u} + (u~u)V t A[\u + (l-\)u}^ +t j d\K(Xt,Xx) 

(2.64) 



A[u] = A(t,x,u,u t ,u x ), 



(2.65) 
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T> t A[u] = A t (t,x,u,u t ,u x ) + u t A u (t,x,u,u t ,u x ) + u tx A Ux (t, x, u, u t , u x ) 

+c(u)(c(u)u x ) x A ut (t, x, u, u t ,u x ), (2.66) 
K(t, x) = {u tt - c(u)(c(u)u x ) x )A[u], (2.67) 

with u = u(t,x) being any function chosen such that the expressions A[u] and u tt — 
c(u)(c(u)u x ) x are non-singular. 

Thus we see that all nontrivial first-order conservation laws of the wave equation 
( [2.47 ) are determined by multipliers of first-order which are obtained as solutions of the 
augmented system of symmetry determining equations ( 2.6C ) and ( [2.61 ). 

We now use the determining system ( 2.60[) and ( 2.61 ) to completely classify all first- 
order conservation laws in terms of corresponding multipliers 

A(t,x,u,u t ,u x ). (2.68) 

The classification results are summarized by the following two theorems. 



Theorem 2.2.1 For arbitrary wave speeds c(u), the only multipliers of form ( 2.6t ) ad- 
mitted by the wave equation (2.4' ) are 



A = u t , A 



u x ,A 



tut + xu x 



(2.69) 



These multipliers define symmetries given by time-translation, space-translation, and 



time-space dilation, respectively, which lead to conservation laws for energy (2.54), mo- 



mentum (2.55), and a dilational quantity. 



Theorem 2.2.2 The wave equation ( 2. 4% for non- constant wave speed ciu) admits ad- 
ditional multipliers of form ( 2.61 ) iff c(u) = cq(u — ito)~ 2 in terms of constants cq,uq. 
For these wave speeds the additional admitted multipliers are given by 

A=t 2 u t -t(u-u ), (2.70) 
A = x 2 u x + x(u — uo), (2-71) 
A = tu t — xu x — (u — uo). (2.72) 



The symmetries defined by these additional multipliers correspond to two conformal 
(Mobius) transformations and one scaling transformation on independent variables t, x, 
accompanied by a scaling and shift of u. The conserved densities $* for the three corre- 



sponding conservation laws are obtained by the construction formula (2.64). We choose 
u = const =/= uq to avoid the singularity at u = uq. This leads to K = 0. Using the linear 



dependence of A(t, x, u, u x , u t ) together with the property A = given by Eqs. ( 2.70 ) 
to (2.72), we find that after integration by parts the formula reduces to 



= -tH(A[tt]+A[u] + ((u-«)A,Ju]) s ) +-(u-u)(A t [u]+A t [u]+tHAJti]) 

' (2.73) 



+~u x 2 c{uyA ut [u] 



up to a trivial conserved density D x 0. Since the terms in Eq. (2.73) are non-singular for 
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any constant u, we can now set u — uq without loss of generality. Then, substituting 



Eqs. (2.70) to (2.72) for A, we obtain respectively 



$* = ^t 2 u t 2 - t(u - u )u t + \{u- u f + U 2 u x 2 c 2 {u - u )-\ (2.74) 
$ = x 2 u x u t + x(u — u )u t , 



1 tu 2 - (XU X + U - U )u t + ^-tu x 2 C 2 (u - Uq) 4 , 



1 

2 ■ v "~" ' ~ ' 2 

which represent two conformal quantities and a scaling quantity. 



(2.75) 
(2.76) 



Proof of Theorems 2.2.1 and 2.2.2: We start from the determining equation (2.61) 



and expand it in explicit form using 

(2.77) 
(2.78) 

This yields 

= 2A +A, +A u t + cc'u 2 A -A -A u x + (c 2 A -A )u xx . (2.79) 

U ZUf UUt Z X UtUt xu x uu x x V U±Ut u x u x / xx \ / 



V t A = A t + A M w t + A Ki « fI + A Ut (c 2 u X2; + cc'm x 2 ), 



Since A does not depend on u xx , Eq. ( 2.7£ ) separates into 

= 2A„ + A tut + A uu u t - A xUx - A uu u x 



(2.80) 
(2.81) 



Next we bring in the symmetry determining equation (2.6C). Expanding it similarly, we 



find after use of Eq. (2.80) that its highest derivative terms involve u xx and u tx . Hence, 
since A does not depend on second-derivatives of u, these terms can be separated, leading 
to 

= (A tUt + A uu u t - A XUi - A uu u x )c + (A u u t + A u u x - A)c' + A UtU u x 2 c 2 c', 

(2.82) 



= A tnx + Au t - (A XUt + A uu u x )c 2 + A v 



To proceed, we combine Eqs. ( p. 81 ) and (2.82) to yield 

°=( A « t ^+ A „^x-A)c'-2A ti c 



(2.83) 



(2.84) 



which is a first order linear PDE for A in u, u t , u x . The solution of Eq. (2.84) is given by 
A = c- 1/2 f{t,x,a,f3), a = c 1/2 u t ,f3 = c 1/2 u x . (2.85) 



Substituting Eq. (2.85) into Eq. (2.80), we obtain c 2 f aa — fpp which separates into 

fact = fpfi = (2.86) 

since c(u) ^ const. Hence we have 

/ = a(t, x)a + b(t, x)(i + g(t, x)a(i + h(t, x). (2.87) 
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Next, we substitute Eqs. ( |2.87| ) and ( |2.85| ) into Eq. ( |2.83| ) to obtain 



= 6 t - c 2 a x + -c -1 /V(t*t 2 + c 2 U:E 2 ).g. 



This immediately yields 
Thus, so far we have 



g = 0, b t = a x = 0. 
A = a(t)u t + b(x)u x + h(t, x)c(u)~ 1/2 . 



Now, using Eq. ( |2.90| ) we see that Eqs. ( p.81| ) and fl2.82|) both reduce to 

= -cr 3/2 ch + a t -b x . 



(2.88) 
(2.89) 
(2.90) 

(2.91) 



This leads to two cases to consider. If the wave speed c{u) is arbitrary, then we must 
have 



Case (i): 



a t = b x , h = 0. 



(2.92) 



Alternatively, the only other possibility is for the wave speed c(u) to satisfy a first-order 
ODE, so we then have 



Case (ii): 



,-3/2 / = - 



c = c = const, ch = a t — b x 



(2.93) 



Note in this case we require c ^ in order for c(u) to be non-constant. 

Finally, we return to the symmetry determining equation ( [2.60 ) and consider the terms 
that remain after we substitute Eq. (2.90). The analysis proceeds according to the two 
cases. 

In case (i), the remaining terms in Eq. (2.60) reduce to 

= -c 2 b xx u x + a tt u t 

which separates into 

b xx — a tt = 0. 

Solving this equation and using Eqs. (HH) and ( ^89|) , we obtain 

a = ao + ait, b = bo + a±x 



where ao,6oi^i are constants. Consequently, from Eq. (2.90), we have 

A = agUt + bou x + a\{tu t + xu x ). 



(2.94) 
(2.95) 
(2.96) 

(2.97) 



Hence, these are the only multipliers of first-order admitted by the wave equation ( 2.47 ) 
for arbitrary wave speeds c(u). This establishes Theorem 



2.2.1 



Finally, in case (ii), using Eq. (2.95) for h and Eq. (2. 



for a and 6, then eliminating 



d by Eq. ( [2.93 ), we find that the remaining terms in Eq. ( 2.60 ) simplify considerably to 

= c~ 1 / 2 h tt -c~ 3/2 h xx . (2.98) 
Hence, since c(u) is non-constant, we obtain 

htt = h xx = 0. (2.99) 
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Now we solve for a, b by combining Eq. ( |2.99 ) with Eqs. (2.93) and (2.89), which gives 

am = b xxx = 0. (2.100) 

Hence we have 

a = ao + ait + a 2 t 2 , b — bo + b\x + b 2 x 2 , (2.101) 

(2.102) 



where ao, ai, 02, &0: 01, 62 are constants. Then Eq. ( 2.93 ) yields 

ch = a\ + 2a2^ — 61 — 262X. 



Consequently, from Eq. ( 2.90 ), we obtain 

A = a u t + b u x + ai(tu t + c -1 / 2 /c) + b\(xu x - c~ 1 ^ 2 /c) 
+a 2 (t 2 u t + 2tc~ 1/2 /c) + b 2 (x 2 u x - 2xc" 1/2 /5) 

where the wave speed is given from Eq. ( 2.93| ) by 

c(u) = cq(u - u )~ 2 

with c = -2co_ 1/2 , c > 0. 

In Eq. (2.103) we have a six parameter family of admitted multipliers. The parameters 
ao,bo,ai = b± yield the three multipliers of case (i). Hence the present case has three 
additional multipliers arising from the parameters ai = —b±, a 2 , b 2 - This establishes The- 
orem [2.2.2j . □ 



(2.103) 



(2.104) 



2.3 Klein-Gordon wave equations 

Consider the class of Klein-Gordon wave equations 

G = ut x - g{u) = 



(2.105) 



with a nonlinear interaction g(u). This class has a variational principle given by the 
action 



S = — I {—u t u x + h(u))dtdx, h'(u) — g{u). 



(2.106) 



Since the general Klein-Gordon equation (2.105) is variational, its symmetries with in- 
finitesimal generator Xit = r) || Q satisfy the determining equation 



= D t D x r] - g'(u)rj when G = 



(2.107) 



which is self-adjoint. Hence the adjoint symmetries of Eq. ( 2.105 ) are the solutions of 

Eq. ( |2.107D . 

The symmetries Xu = 77 admitted by the general Klein-Gordon equation (2.105) clearly 
consist of t and x translations 

7]=U t ,T)=U x , (2.108) 

as well &s &t,x boost 

7] = tu t — xu x . (2.109) 
These are easily checked to be symmetries of the action, leaving S invariant up to a 
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boundary-term, with 



XS 



D t (—u t u x + h(u))dtdx, XS 



D x (-u t u x + h{u))dtdx 



(2.110) 



under the respective translations Xu = u t and X« = u x , and with 



XS 



^ — D t (—tu t u x + th(u)) + D x (—xutu x + xh{u))*jdtdx (2.111) 



under the boost Xu = tu t — xu x . By Noether's theorem, combining the invariance of the 
action and the general variational identity 



XS 



([u tx - g{u))rj - D t (^rju x ) - D x (-r]Utjjdtdx, Xu = rj, 



(2.112) 



we obtain corresponding local conservation laws D t $ +D x § x = on all solutions u(t, x) 
of the Klein- Gordon equation ( [2.105 ). The conserved densities are given by 

& = h(u),$ x = -^u t 2 , (2.113) 
1 



-u x \<S> x = h(u) 



from the translations, and 



1 



$* = -xu x 2 +th(u),<i> x 



1 



xh(u), 



(2.114) 



(2.115) 



2 x " 2 

from the boost. These comprise all of the first-order local conservation laws for the general 
Klein- Gordon equation ( 2.105 ). 



The Klein-Gordon equations in the class (2.105) include 

u tx = sin it, sine-Gordon equation 

utx = e" ± e _tl , sinh-Gordon equation 

Utx — e u ± e~ 2 ", Tzetzeica equation 

utx = e 11 - Liouvillc equation 



(2.116) 
(2.117) 
(2.118) 
(2.119) 



The first three are soliton equations while the last is a linearizable equation. These are 
singled out [Q as nonlinear wave equations that are known to be integrable in the sense 
of admitting an infinite number of higher-order local conservation laws |7], [l0| . For each 
equation the conservation laws fall into two sequences where $' and ® x depend purely 
on u and t derivatives of u in one sequence, and purely on u and x derivatives of u in 
the other sequence (corresponding to the reflection symmetry t x). 

Here, for the class of nonlinear Klein-Gordon equations (2.105), we consider local 
conservation laws of higher-order 

-t 







(2.120) 



with $ , <& x depending either purely on t, u, and t derivatives of u, or purely on x, u, and 
x derivatives of u, on all solutions u(t, x) of Eq. ( [2.105 ). 

All nontrivial conservation laws ( |2.120 ) with conserved densities of the particular form 



$ (x, u, d x u, . . . , d%u), & x (x, u, d x u, . . . , d%u) 



(2.121) 
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can be constructed from multipliers A on the Klein-Gordon equation ( 2.105| ), analogous 
to integrating factors, with the dependence 

K(x,u,d x u,...,8Pu) (2.122) 

where p — 2q — 1. In particular, moving off the Klein-Gordon solution space, we have 

D t & + D x $ x = (u tx - g(u))A Q + D x (u tx - 5 (u))A 1 + • • • (2.123) 

for some expressions A Q , A 1; . . . with no dependence on u t and differential consequences. 
This yields (after integration by parts) the multiplier 

D^ + DJO? -T) = (u tx -g{u))k, A = A - D x A x + • • ■ (2.124) 

where r = when u is restricted to be a Klein-Gordon solution. There are corresponding 
conservation laws of mirror form produced by the transformation x <-» t, d x u <-> d^u, with 
<F <-> $ x . We now derive an augmented symmetry determining system which completely 



characterizes the conservation law multipliers ( 2.122 ). 

The determining condition for multipliers A is that (u tx — g(u))A must be a divergence 
expression for all functions u(t,x) (not just Klein-Gordon solutions). For multipliers of 
the form ( 2.122) ), this condition is expressed by 



0= E u ((u tx -g)h) 

= D t D x A - g'A + A u (u tx - g) - D x (A dxU (u tx -g)) + ...+ {-D x f{A dlu (u tx - g)) 

(2.125) 

where E u = d u — D f d u — D,B U + D f DB u + D^du + D 2 t d u + ■ ■ ■ is the stan- 

1 t X l J- i x X xx l f-f- ^^^^^^ 

dard Euler operator which annihilates divergence expressions. Equation ( 2.125 ) is lin- 
ear in u tXl Ut xx , . • ., and hence splits into separate equations. We organize the splitting 
in terms of G,D x G, . . . ,D p x G. This yields a split system of determining equations for 

A(x, u, d x u, . . . , d^u), which is found to consist of the Klein-Gordon symmetry determin- 
ing equation on A 

= V t D x A - g'(u)A (2.126) 
and extra determining equations on A 

° = 2A„ + E(-^) feA ^' ( 2 - 127 ) 

k=2 

o = (i + (~m\ iu + ((-iy-j-i)D x A d ^ u+ ± 1 u^-y(-D x ) k - j A dkxU , 

k=j+2 J ' y J '' 

i = l,...,p-l (2.128) 
= (l + (-l)")A flSlt . (2.129) 

Here T> t = dt + u t d u + g(u)d Ua . + g'(u)u x d Ux:l . + ■ ■ ■ is the total derivative operator which 
expresses t derivatives of u through the Klein-Gordon equation ( [2.105 ). Consequently, 



one is able to work on the space of Klein-Gordon solutions u{t, x) in order to solve the 



determining system (2.126) to (2.129) to find A(x,u,d x u, . . . ,d^u) 
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The solutions of the determining system ( |2.126 ) to (2.129) are the multipliers (2.122) 
that yield all nontrivial conservation laws of the form (2.121). Noether's theorem estab- 
lishes that there is a one-to-one correspondence between symmetries of the action and 
multipliers for nontrivial conservation laws of the Klein-Gordon equation ( |2.105 ). Conse- 
quently, it follows that the extra determining equations (2.127) to (2.12E) represent the 



necessary and sufficient condition for a symmetry to leave the action ( 2.106 ) invariant 
up to a boundary term. 

The explicit relation between multipliers A and conserved densities $ ,<S> X for conser- 



vation laws (2.121) is summarized as follows. Given a conserved density $ , we find that 
a direct calculation of the utx terms in D t <& + D x <& x yields the multiplier equation 



D t & + D x & = (u tx - g(u))E Ua: ($') + D X T 



(2.130) 



where E Ux — d u — D x d u + • • • is a restricted Euler operator, and T is proportional to 
Utx ~ g{u) and differential consequences. Thus we obtain the multiplier 

A=K x ($ t ). (2.131) 



Conversely, given a multiplier A, we can invert the relation (2.131) using Eq. (2.130) to 
obtain the conserved density 



d\(u x — u x )A[Xu + (1 — X)u] 



(2.132) 



where A [it] = A(a;, u, d x u, . . . , d^u), and u — u(x) is any function chosen so that the ex- 
pressions A [it] and g(u) are non-singular. In particular, if A[0] and g(0) are non-singular, 
then we can choose u = 0, which simplifies the integral ( [2.132 ). 

Thus, from Eqs. ( 2.131 ) and ( 2.132 ), we see that all nontrivial Klein-Gordon conserva- 
tion laws ( 2.121 ) of order q are determined by multipliers ( [2.122 ) of order p = 2q—l which 
are obtained as solutions of the augmented system of symmetry determining equations 
( [2.126D to Q2.129| ). 

Through the determining system fl2.126| ) to ([2.129D , we now derive a classification 
of nonlinear Klein-Gordon interactions g(u) that admit higher-order conservation laws 
( [2.121 ) starting at order q = 2. The classification results are summarized by three theo- 
rems. 



Theorem 2.3.1 The nonlinear Klein-Gordon equation (2.106) admits a conservation 



law (2.121) of order q — 2 iff the nonlinear interaction is given by one of 

g(u) = e a , g(u) = sinw, g(u) = e u ± e~ u (2.133) 
to within scaling oft,x,u, and translation of u. 

(If we allow u to undergo complex-valued scalings and translations then the Klein-Gordon 
equations for g{u) = sinu and g{u) — e" ± e~ u are equivalent.) 



The Klein-Gordon equations arising from Theorem 2.3.1 are the Liouvillc equation 



( [2.119 ), sine-Gordon equation (2.116) and related sinh-Gordon equation ( 2.117| ). (The 
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Tzetzeica equation ( |2.118 ) is absent in this classification because its first admitted higher- 
order conservation law (2.121) is of order q — 3.) Since each of these Klein- Gordon 
equations is known to admit an infinite sequence of higher order conservation laws (2.121) 
for q > 2, this leads to an integrability classification. 



Corollary 2.3.2 A nonlinear Klein-Gordon equation ( 2.1 PL ) is integrable in the sense 
of having an infinite number of conservation laws (2.121) up to arbitrarily high orders 
q > 1 if it admits one of order q = 2. 



by the Klein- Gordon equation ( 2.10\) with nonlinear interactions ( 2.13l 



Theorem 2.3.3 The multipliers for the second-order conservation laws (2.121 ) admitted 

are given by 

(2.134) 
(2.135) 



A U xxx -f- ^ ^X ? 

1 3 

U XXX ~ U x , 



for the sine-Gordon equation (2. lit ) and sinh-Gordon equation (2.117), respectively, and 
A = (u xxx - u x u xx )f(: + f x + u x f = D x f + u x f, (2.136) 
depending on an arbitrary function 



f = /(&,£)> t = Uxx- ' 



(2.137) 



for the Liouville equation (2. Hi ). 



From the construction formula ( p. 132 ) the corresponding conserved densities for the 
multipliers (2.134) and ( 2.135] ), respectively, are given by 



zb —U, 



4 + dj, 



(2.138) 



where we have used u = 0, since A[0] =0 and g(0) = const. Similarly, the conserved 
density for the multiplier ( 2.136| ) is given by 

$ = (u x - u x )D x \ d\f(x, \{u xx - u xx ) ~ -\ 2 {u x - u x ) 2 ) 



+ (u x 



f / d\ Xf(x, X(u xx - u xx ) - -X 2 (u x - u x ) 2 ) 



f(x,Z)d£ + D x t 



(2.139) 



where now u = u{x) is any function chosen so that A[u] = D f(x,£) + u x f(x,^) and 
g(u) = exp(u) are non-singular expressions. 

We remark that the existence of the conservation law (2.139) involving an arbitrary 
function reflects the well-known classical integrability || (in the sense of explicit inte- 
gration) of the Liouville equation. 



Proof of Theorem 



2.3.1 
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From Eqs. (2.121) and ( [2.122| ), conservation laws with u, u x , u xx ) and 
<& x (x, u, u x , Uxx) of order q = 2 correspond to multipliers 

A(x,u,u x ,u xx ,u xxx ). (2.140) 
We start with the extra determining equations (2.127) to ( 2.129| ). These become 

° = 2A U +£X^-^A W (2.141) 

= D x K xx -D 2 x A Uxxx , (2.142) 

" V, "A, . • (2.143) 



Note that Eq. (|2.142|) is a differential consequence of Eq. ( |2.143| ), while Eq. fl2.14l|) 
combined with a differential consequence of Eq. ( [2.142 ) reduces to 



= A, 



Consider Eq. (EEl43~|). Its highest-order term is A„ „ u xxxx , and hence A„ 
This yields, using Eq. ( gig ), 

A U x j ^xx^^xxx H~ b(^X, Tlxi ^xx) ■ 



Then the remaining terms in Eq. (2.143) give 



^U xx — ^0? ~T~ ^Ux^XX' 



(2.144) 

tra = 0. 

(2.145) 
(2.146) 



We now turn to the symmetry determining equation ( 2.126| ). This becomes, taking 
into account Eqs. fl2.140p and ( |2. 144f ) , 



= V t A x + u xx V t A u + u xxx V f A u 



^X X X x ^ t ^ u x 



-KD x g + A u D 2 x g + AD 3 x g - A,/. 



(2.147) 



We substitute Eq. (2.145) into Eq. (2.147) and separate it into highest derivative terms 
in descending order. The terms of highest order involve u xxxx , which yield 



= au*g + a,u xx u x g' . 



(2.148) 



Since g depends only on u, while a does not depend on u, the terms in Eq. (2.148) 
can balance only in two ways: either g and g 1 are proportional, or g and g' are linearly 
independent with vanishing coefficients. This leads to two cases to consider. 

Case (i): g,g' linearly independent 

In this case a Ux — a Uxx = 0, so a = a(x). Consequently from Eq. ( |2.146 ), we have 
b Uasx = a Xl and hence 

b = a'(x)u xx + c(x,u x ). (2.149) 

To proceed, we return to the symmetry determining equation ( BT47 ). The highest order 
remaining terms now involve u xx , which lead to 

= a g + 3au x g" + c UxUx <7. (2.150) 

By taking a partial derivative with respect to u x we obtain = Sag" + c UxUxUx g. Since 
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g depends only on u, while a does not depend on u, this equation immediately yields 

g" = ag, a = const, (2.151) 
Cu,«,u, = Sao, (2.152) 



(2.153) 



where a ^ in order for g(u) to be nonlinear. From Eq. (2.152) we obtain 
c = --aa(x)u x 3 + c 2 (x)u x 2 + ci{x)u x + c a (x). 

Now we find that Eq. ( 2.150| ) reduces to = a'g' + 2c 2 g, which immediately separates 
into the equations 

a' = c 2 = 0. (2.154) 



Finally, the terms remaining in the symmetry determining equation (2.147) simplify to 
= c\ x g — cog' , and thus we have 



co = CXx = 0. 

From Eq. fl2.145| ), Eq. Q2.149D , and Eqs. ( [2.153D to @.155] ), we obtain 



A = a(u x 



— au x ) + ciu x , a — const, c\ — const. 



(2.155) 
(2.156) 



This multiplier is admitted for all nonlinear interactions g(u) satisfying Eq. (2.151) with 
g(u) and g'(u) being linearly independent. The general solution for g(u) breaks into two 
forms: 



g(u) = a sin(^/\a\u + /3), if a < 0; 
g(u) = a{e^ u+0 ± e"^ 11 ^), if a > 0; 



(2.157) 
(2.158) 



with arbitrary constants a, f3, o ^ 0. By a scaling and a translation -^/jcrju + f3 — > u, and 
a scaling a-^/jo^fi — > t, we see that the resulting Klein-Gordon equation ( |2.105 ) becomes, 
respectively, the sine-Gordon equation ( 2.116| ) or the sinh-Gordon equation (2.117). This 
completes the classification in case (i). 



Case (ii): g,g' proportional 
In this case we have 



fJ 



ag,a — const 



where a ^ in order for g(u) to be nonlinear. Now Eq. (2.148) becomes 

= a Ux + aa Uiic u x . 
This is a linear first-order PDE for o in m x , u xx . The solution is 

a = a(x, f), £ = u xx ~ \^ju x 2 . 
Then from Eq. ( |2.146| ) we obtain 

b = a x (x,£) + au x (a(x, £) - a(x,£)u xx ) + c(x,u x ), a= ad£,. 



(2.159) 
(2.160) 

(2.161) 
(2.162) 



To proceed, we return to the symmetry determining equation (2.147). We find that 
the highest order terms involving u xxx 2 and u xxx vanish as a consequence of Eq. (2.160). 
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Then we find that the next highest order terms which involve u xx in Eq. (2.147) yield 
Cua-Ua, = 0. Hence, we obtain 

c = c\{x)u x + cq(x). (2.163) 
Finally, the remaining terms in Eq. ( [2.147 ) reduce to 

= c[-acQ. (2.164) 

Thus, from Eq. fl2.145| ), Eqs. fl2.16l| ) and ( |2.162| ), Eqs. ( [2.163D and ( |2.164D , we obtain 



A = a^u xxx + <Ju x (a - u xx a^) + a x + c\u x -\ c[, £ 



(2.165) 



with arbitrary functions c\{x), a{x,^). This multiplier is admitted for all nonlinear inter- 



actions g(u) satisfying Eq. (2.159). The general solution of this equation is 



(2.166) 



with arbitrary constants j3,<j =/= 0. By a scaling and a translation au + (3 — > u, and a 



scaling at — » t, we see that the resulting Klein-Gordon equation ( 2.105 ) becomes the 
Liouville equation. This completes the classification in case (ii). □ 
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